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Abstract 

After the recent approval by the Italian Space Agency (ASI) of 
the LARES mission, which will be launched at the end of 2008 by 
a VEGA rocket to measure the general relativistic gravitomagnetic 
Lense- Lhirring effect by combining LARES data with those of the ex- 
isting LAGEOS and LAGEOS II satellites, it is of the utmost impor- 
tance to assess if the claimed accuracy < 1% will be realistically ob- 
tainable. A major source of systematic error is the mismodelling 5 Jg in 
the static part of the even zonal harmonic coefficients Ji,£ — 2,4,6, .. 
of the multipolar expansion of the classical part of the terrestrial grav- 
itational potential; such a bias crucially depends on the orbital config- 
uration of LARES. If for S Jg the difference between the best estimates 
of different Earth's gravity solutions from the dedicated GRACE mis- 
sion is conservatively taken instead of optimistically considering the 
statistical covariance sigmas of each model separately, as done so far 
in literature, it turns out that, since LARES will be likely launched in 
a low-orbit (semimajor axis a < 7600 km), the bias due to the geopo- 
tential may be up to ten times larger than what claimed, according to 
a calculation up to degree £ = 20. Taking into account also the even 
zonal harmonics with £ > 20, as required by the relatively low altitude 
of LARES, may further degrade the total accuracy. Should a nearly 
polar configuration (inclination to the Earth's equator z w 90 deg) be 
finally implemented, also other perturbations would come into play, 
further corrupting the measurement of the Lense-Thirring effect. The 
orbital configuration of LARES may also have some consequences in 
terms of non-gravitational perturbations and measurement errors. 
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1 Introduction 



The Italian Space Agency (ASI) recently made the following official an- 
nouncement pQ: "On February 8, the ASI board approved funding for the 
LARES mission, that will be launched with VEGAs maiden flight before 
the end of 2008. LARES is a passive satellite with laser mirrors, and will 
be used to measure the Lense-Thirring effect." The italian version of the 
announcement yields some more information specifying that LARES, de- 
signed in collaboration with National Institute of Nuclear Physics (INFN) , 
is currently under construction by Carlo Gavazzi Space SpA; its Principal 
Investigator (PI) is I. Ciufolini and its scientific goal is to measure at a 
1% level the general relativistic gravitomagnetic Lense-Thirring effect [2], 
known also as frame-dragging, in the gravitational field of the Earth. 

After what recently happened with the much more expensive Gravity 
Probe B (GP-B) [3j[I], which should have measured another gravitomagnetic 
effect in the terrestrial gravitational field, i.e. the Schiff precession [5j of four 
gyroscopes carried onboard at an accuracy much worse than the expected^ 
1%, the success of the cheaper LARES mission would be of great scientific 
significance. 

The major source of systematic uncertainty in such kind of satellite- 
based measurements is represented by the impact of the static part of the 
even (I = 2,4,6...) zonal (m = 0) harmonic coefficients Ji of the multipolar 
expansion of the classical part of the terrestrial gravitational potential [7] 
which induce noising effects having the same signature of the relativistic 
secular precessions of the satellite's node £1 and a much larger amplitude. 
Minimizing such a bias is of crucial importance and it strongly depends on 
the orbital parameters of LARES, especially the semimajor axis a and the 
inclination i to the Earth's equator. When it was proposed for the first time 
with the name of LAGEOS III [3|9] the presently known LARES [TO] had the 
same semimajor axis of the existing LAGEOS satellite (a = 12270 km) and 
inclination supplementary to LAGEOS («lageos = 110 deg, Glares = 70 
deg) ; such a configuration would have allowed to use the sum of the nodes 
of both satellites enhancing the Lense-Thirring signal and cancelling out, at 
least in principle, the competing effect of all the even zonals. The advent of 
recent models of the Earth's gravity field by the dedicated missions CHAMP 
[TT] and GRACE [12J and the idea of combining LARES data with those 
of LAGEOS and LAGEOS II [TOJ Q3j, according to an idea put forth in 
|17j . suggested to abandon the original stringent requirements allowing for 

^ee 0. 
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a lower orbit and a less narrow range for the inclination |14j . 

No details at all are released concerning the orbit in which LARES will 
be finally injected: ASI website says about VEGA that [T5]: "[...] VEGA 
can place a 15.000 kg satellite on a low polar orbit, 700 km from the Earth. 
By lowering the orbit inclination it can launch heavier pay loads, whereas 
diminishing the payload mass it can achieve higher orbits. [...]" In the 
latest communication to INFN, Rome, 30 January 2008, [16] I. Ciufolini 
writes that LARES will be launched with a semimajor axis of approximately 
7600 km and an inclination between 60 and 80 deg. 

As it will be shown, a detailed and realistic evaluation of the impact of 
the bias due to the geopotential shows that such a low-orbit configuration 
should not allow to reduce such a systematic noise below the 1% level being, 
instead, a few times larger. 



2 An evaluation of the systematic bias due to the 
even zonal harmonics of the geopotential 

The Lense-Thirring effect consists of a tiny secular precession of the node f2 
of the orbit of a satellite moving around a central slowly rotating mass 

■ _ 2GS 

LT " c 2 a 3 (1 _ e 2 )3 /2' W 

where G is the Newtonian constant of gravitation, S is the spin angular 
momentum of the central body, a and e are the semimajor axis and the 
eccentricity, respectively, of the satellite orbit. For the LAGEOS satellites 
it amounts to about 30 milliarcseconds per year (mas yr _1 ). The observable 
to be used to measure it is the following linear combination of the nodes' 
rates [H] O of LAGEOS, LAGEOS II and LARES 

^LAGEOS + ^LAGEOS II + ^LARES. (2) 

5fl is an Observed-minus-Calculated (O-C) quantity for the satellites' nodal 
rates to be observationally constructed by processing the laser-ranging data 
without modelling the gravitomagnetic force. It accounts, among other 
things, for all the unmodelled (like the Lense-Thirring effect itself) or mis- 
modelled dynamical effects affecting the spacecraft nodes being equal to 
Alt + A) where A represents all the classical mismodelled/unmodelled ef- 
fects of gravitational and non-gravitational origin. The combined Lense- 
Thirring signature, which is expected to be present in the signal of eq. (|2|), 
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IS 

X LT = ^MiEOS + c ^LAGEOS II + ^fiLAHES (3) 
The coefficients c\ and c 2 entering eq. (|2|) and eq. (J3j) ar e denned as 



QLARESQLAGEOS_f2LAGEOSQLARES 

ci — 



qLAGEOS IIqKARES_j-2LARESqLAGEOS II • 



qLAGEOSqLAGEOS I^j^LAGEOS IIqLAGEOS 
C 2 qLAGEOS IlQLARES_f : 2LARESQLAGEOS II 1 



(4) 



The quantities entering eq. @ are defined, for the satellite x, as 

qx _ Ui 'class /r\ 

where 

^class = (6) 
1=2 

is the classical node precession [7] induced by the even {I = 2, 4, 6, ...) zonal 
(m = 0) harmonic coefficient^! Jg of the multipolar expansion of the New- 
tonian part of the terrestrial gravitational potential; for example, for I = 2 
we have 

3 / R\ cos i 
2 U a fW 



where R is the equatorial radius of the central body of mass M and n = 
\J GM/a 3 is the satellite's Keplerian mean motion. The coefficients (Vg, and 
thus also ci and C2, are functions of the semimajor axis a, the eccentricity 
e and the inclination i of the satellite x. The combination of eq. (|2|), along 
with eq. (H|), is designed, by construction, to cancel out the mismodelled 
parts 5J2 and 5J4 of the static and time-varying components of the first two 
even zonals, being affected by all the other ones of higher degree 5Jq,5J$, ■■■ 
The evaluation of the systematic percent error 5u due to the geopotential 
has been so far performed [131 OS E! according tqj 

/ EZ2 I ^ AGEOS + ci^ AGEOS 11 + c 2 ^ ARES I 5Ji \ 
6fl - I f^GEOS + ^LAGEOS II + ^LARES J W ° (8) 

They are defined as Je — —\/2£ + 1 Cto, where C10 are the normalized gravity coefH- 
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cients 

3 



Since we are dealing with a test of fundamental physics, one has to be quite conser- 
vative. This is the reason why we do not consider a Root-Sum-Square evaluation of the 
systematic error and prefer to linearly adding the individual terms. 
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Table 1: Systematic percent error 5/i in the measurement of the Lense- 
Thirring effect with the combination of eq. ([2]) according to eq. ([8]) and 
the calibrated covariance sigmas u^? of the Earth's gravity model EIGEN- 
GRACE02S up to degree I = 20 for different values of the inclination i of 
LARES (q = 7600 km, e = 0.001). 





i = 60 deg 


% = 70 deg 


i = 80 deg 


5n (EIGEN-GRACE02S) 


4.3% 


1.6% 


3.8% 



by taking the covariance sigmas a-Q m of a given Earth's gravity solution and 
using them for the uncertainty 5Ji in the even zonals; in this way it has 
been always said that model X yields a total error x%, model Y yields a 
total error y% and so on. See, e.g., Table Q] for the EIGEN-GRACE02S 
solution. Unfortunately, it is difficult to assess how realistically the sigmas 
of the covariance matrix of a given solution reflect the true uncertainties in 
the even zonals: in some cases, they are the mere formal, statistical errors, 
in other cases they are calibrated. However, also their calibration is often 
difficult to be realistically implemented, coming from a rather ad-hoc proce- 
dure. It is important to understand what kind of covariance information is 
taken for a particular gravity field solution. A filter covariance, i.e. the un- 
certainties are given as they are determined by the estimation filter, tends 
to yield rather optimistic errors. A more reliable approach— not followed 
so far in the publicly available models due to its great difficulties— would 
consist in using the so-called consider covariance matrix, i.e. uncertainties 
of modelling parameters like, e.g., the Earth's rotation, are taken into ac- 
count. The consider covariance would give a rather good indication of the 
real uncertainties. 

A much more realistic and quantitative approach, which allows every- 
one to make an own idea, consists, instead, in taking for bJg the differences 
ACffl = \Ceo(X) — Ceo(Y)\ between the best estimates of the model X and 
the model Y, which in many cases, are significantly larger than the sum of 
the sigmas o-q (X) + cr^? (Y). Such an approach is usually applied by re- 
searchers involved in gravity field determination [19] and has been followed 
in the case of the previous tests of the Lense-Thirring effect with the LA- 
GEOS satellites by Ciufolini himself [201 QU by comparing the GEMT-3S 
and JGM-3 gravity fields. See also [21 1 . In this paper we will follow it for 
several published models, retrievable al|j [22] along with references, obtained 

4 Concerning the model GGM03S [53], I gratefully thank J. Ries, CSR, for having 
provided me with its spherical harmonic coefficients. 
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by different institutions from GRACE data. Concerning the choice of the 
models, we did not include in our analysis the latest solutions by GFZ en- 
compassing data from LAGEOS as well because of the strong possibility 
of a-priori 'imprint' of the gravitomagnetic signature itself. Moreover, in 
view of the fact that the latest LAGEOS satellites SLR analyses concern- 
ing the Lense-Thirring effect measurement (but not only in this case) have 
been based on satellite-only data, without considering the contribution of 
gravimetry and altimetry surface measurements, we did not considered so- 
lutions including such kinds of data like, e.g., EIGEN-CG03C. It has been 
decided to use for LARES a = 7600 km and a = 7000 km, and three differ- 
ent values for the inclination: i = 60, 70, 80 deg. As can be noted, claiming 
a bias < 1% is unrealistic. The results are graphically depicted in Figure Q] 
(a = 7600 km) and Figure M (a = 7000 km); see also Tabled 

Analogous analyses conducted for a = 7600 km and a = 7000 and i = 85 
deg, i = 89 deg, i = 91 deg, which corresponds to the most likely orbital 
configuration allowed by VEGA, show a neat worsening, as shown by Table 
® 

In Figure [3] and [H we depict the situation for a = 7600 — 7000 km, 
e = 0.001, 85 deg < i < 89 deg. 

A caveat concerning the results presented so far is that, in principle, they 
might be optimistic because they have been obtained by only accounting for 
the even zonals up to degree t = 20. But, with such a small semimajor 
axis, it maybe required to take into account also the impact of the other 
even zonals of higher degree I 20. Such doubts are, in fact, enforced by 
a straightforward calculation according to the standard approach by [7] for 
a = 7600, km i = 71 deg, e = 0.001 whose results are presented in Tabled! 

The only valid solution would be to use the originally proposed orbital 
configuration for LARES (a = 12270 km, i = 70 deg), as shown by Figure 
[5j as can be noted, for a deviation of 1 deg from the optimal choice i = 70 
deg the systematic error would be well below 1%. 

3 The coefficients of the combination 

The inclination in which LARES will be finally inserted into its orbit is 
also crucial in determining the size of the coefficient C2 with which its node 
enters the combination of eq. ([2]). It does matter because it may enhance 
all the perturbations not cancelled out by the coefficients of eq. @). They 
are the non-gravitational ones, whose impact on the node of LARES should 
be reduced to ~ 0.1% of the Lense-Thirring effect thanks to its careful 
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Systematic bias due to the even zonals according to several Earth gravity solutions: a = 7600 km, e = 0.001 




Figure 1: Systematic percent error 5(i due to the even zonal harmonics up 
to degree i = 20 for a = 7600 km, e = 0.001, 60 deg < i < 80 deg according 
to several Earth's gravity models. The uncertainty 5Je has been evaluated 
by taking the differences of the best estimates of the even zonals for several 
pairs of models. 
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Systematic bias due to the even zonals according to several Earth gravity solutions: a = 7000 km, e = 0.001 
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Figure 2: Systematic percent error b\i due to the even zonal harmonics up 
to degree I = 20 for a = 7000 km, e = 0.001, 60 deg < i < 80 deg according 
to several Earth's gravity models. The uncertainty 5Je has been evaluated 
by taking the differences of the best estimates of the even zonals for various 
pairs of models. 
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Systematic bias due to the even zonals according to several Earth gravity solutions: a = 7600 km, e = 0.001 
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Figure 3: Systematic percent error 5(i due to the even zonal harmonics up 
to degree £ = 20 for a = 7600 km, e = 0.001, 85 deg < i < 89 deg according 
to several Earth's gravity models. The uncertainty 5Je has been evaluated 
by taking the differences of the best estimates of the even zonals for various 
pairs of models. 
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Systematic bias due to the even zonals according to several Earth gravity solutions: a - 7000 km, e = 0.001 
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Figure 4: Systematic percent error 5(i due to the even zonal harmonics up 
to degree I = 20 for a = 7000 km, e = 0.001, 85 deg < i < 89 deg according 
to several Earth's gravity models. The uncertainty 5Je has been evaluated 
by taking the differences of the best estimates of the even zonals for various 
pairs of models. 
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Figure 5: Systematic percent error 5(i due to the even zonal harmonics up 
to degree i = 20 for a = 12270 km, e = 0.001, 68 deg < i < 72 deg according 
to several Earth's gravity models. The uncertainty 5Je has been evaluated 
by taking the differences of the best estimates of the even zonals for various 
pairs of models. 
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construction, and the i = 2, m = 1 component of the solar K\ tide whose 
node harmonic perturbation has the same period of the satellite's node [23] 
and for which, of course, the physical properties of LARES are of no concern. 
It turns out that the case for 60 deg < i < 80 deg does not pose problems 
because c\ < 0.5 and c 2 is of the order of 10~ 2 . 

The situation is quite different for inclinations close to i = 90 deg for 
which C2 diverges. Figure [6] shows that for a = 7000 km and 89.0 deg 
< i < 89.9 deg, i.e. the most likely orbital configuration in which VEGA 
will deploy LARES, the coefficient C2 of the node of LARES becomes as 
large as 5, getting even larger (c2 < 55) for 89.90 deg < i < 89.99 deg. 
For a = 7000 km and 89.0 deg < i < 89.9 deg the period of the node of 
LARES, which is the same of the node harmonic perturbation induced by 
the t = 2, m = 1 constituent of the solar K\ tide, not cancelled out by the 
combination of eq. ([2]), amounts to several years, as shown by Figure [7J this 
fact would induce a serious systematic error because the K\ tide would act 
as a superimposed corrupting linear trend over the time spans of some years 
which would be typically used for the data analysis. 



4 The impact of the inclination errors 

Another source of systematic error which may become important for certain 
orbital configurations is the cross-coupling among J2 and the errors in the 
inclinations of the satellites entering the combination of eq. (|2|). Indeed, 
such a bias <5/Ui nc ii is not cancelled out by the coefficients of eq. ([4]) being 
equal to 



^Mincli < 



—Cl L 2 J2 taniL^L 


+ 


-ciCl L 2 II J 2 tan i L n^L 11 


+ 


-c 2 Si I 2 R J2 tan «LR^LR 




OLAGEOS + ClS 7LAGEOS II + C2 QLARES 





(9) 

If we assume 5i ~ 5r/a, with 5r ~ 1 cm for all the satellites, it turns out 
that eq. ([9|) yields a further bias which is about 1% for 60 deg < i < 80, 
but it may become larger for departures of 1 deg from i = 90 deg, as shown 
by Figured) 
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Coefficients c, of LAGEOS II and c„ of LARES 



- LAGEOS II 
LARES 



/(deg) 



Figure 6: Coefficients c\ and C2 of LAGEOS II and LARES, according to 
eq. ©, for a = 7000 km, e = 0.001 and 89.0 deg < i < 89.9 deg. For 89.90 
deg < i < 89.99 deg C2 gets as large as 55. 
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Node period of LARES 



I I I I I I I I I I I I 

89.05 89.1 89.15 89.2 89.25 89.3 89.35 89.4 89.45 89.5 89.55 89.6 89.S5 89.7 89.75 

f(deg) 



89.8 89.85 89.9 



Figure 7: Period of the node of LARES, in years, for a = 7000 km, e = 0.001 
and 89.0 deg < % < 89.9 deg. 
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Systematic error due to the cross-coupling between J 2 and the errors 5 / in the inclinations 




' 89 89.1 89.2 89.3 89.4 89.5 89.6 89.7 89.8 89.9 90 90.1 90.2 90.3 90.4 90.5 90.6 90.7 90.8 90.9 91 

/(deg) 



Figure 8: Systematic percent error <5/ii nc ii due to the cross-coupling between 
J2 and the errors 5i in the inclinations of the satellites entering the combi- 
nation of eq. $2$). We assumed for LARES a = 7600 km, e = 0.001, 89 deg 
< i < 91 deg. 
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5 Some considerations on the non-gravitational per- 
turbations and the measurement errors 



It is worthwhile noting that also the impact of the subtle non-gravitational 
perturbations will be different with respect to the original proposal because 
LARES will fly in a different and lower orbit and its thermal behavior will 
probably be different with respect to LAGEOS and LAGEOS II. A detailed 
treatment of this important subject is outside the scopes of this paper; 
we will only give some insights calling for deeper analyses by more expert 
researchers in the field. The reduction of the impact of the thermal accel- 
erations, like the Yarkowsky-Schach effects, should have been reached with 
two concentric spheres. However, as explained in [25], this solution will 
increase the floating potential of LARES because of the much higher elec- 
trical resistivity and, thus, the perturbative effects produced by the charged 
particle drag. Moreover, drag will increase also because of the lower orbit 
of the satellite, both in its neutral and charged components; preliminary 
calculation point towards a secular decrease of the inclination of LARES 
which maps into a further, uncancelled bias in the node precessions due to 
J 2 which may be relevant over the typical timescales of the test |26| . Also 
the Earth's albedo, with its anisotropic components, should have a major 
effect. 

Another point which must be considered is the realistic orbit accuracy 
obtainable for LARES. Indeed, at a lower orbit the normal points RMS will 
be probably higher with respect to the present RMS obtained for the two 
LAGEOS satellites (a few mm), as we presently know for the Stella and 
Starlette normal points. Of course, such an accuracy is a function of several 
aspects. 

6 Conclusions 

In this paper we have calculated the systematic error due to the even zonal 
harmonics of the geopotential, up to degree £ = 20, on the measurement of 
the Lense-Thirring effect to be performed with the existing LAGEOS and 
LAGEOS II satellites along with the recently approved LARES, which will 
be launched at the end of 2008 by ASI with a VEGA rocket. 

By taking the differences between several Earth's gravity solutions from 
the dedicated GRACE mission instead of optimistically considering the sta- 
tistical covariance sigmas of each model separately, as done so far, we have 
shown that, with the orbital configuration of LARES which should be im- 
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plemented in such a mission (a = 7600 km, 60 cleg < i < 80 deg), the claim 
of reducing the impact of the mismodelling in the even zonals at a < 1% 
level is optimistic because it may be up to ten times larger according to 
calculations up to degree t = 20. Taking into account the even zonals of 
degree i > 20 whose action must be considered because of the relatively low 
altitude of LARES, may further degrade the total accuracy. Should a nearly 
polar, lower-orbit (a = 7000 km, i ~ 90 deg) is implemented, all the figures 
would get even worse: in this case, also other orbital perturbations would 
contribute to corrupt the outcome of the test. The orbital configuration of 
LARES may also have an impact on some non-gravitational perturbations 
and measurement errors. 

Only the originally proposed LARES configuration (a = 12270 km, i = 
70 deg) would damp the systematic geopotential error to a few percent level. 

In conclusion, according to the present-day state-of-the-art in our knowl- 
edge of the terrestrial gravity field, the ongoing LARES mission, as it seems 
it will be implemented, should not be able to allow for a < 1% measurement 
of frame-dragging. Maybe LARES will be more useful for studies of geodesy 
and geophysics. 
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Table 2: Systematic percent error Sfi in the measurement of the Lense-Thirring effect with the combination of 
eq. ([2]) according to eq. (|8|) and the difference AC|o among the best estimates for the even zonal coefficients for 
several Earth's gravity solutions up to degree t = 20 for different values of the inclination i of LARES and for 
a = 7600 km, a = 7000 km (e = 0.001). See also Figure [Hand Figure [j 





i = 60 deg 


i = 70 deg 


i = 80 deg 


a = 7600 km 


a = 7000 km 


a = 7600 km 


a = 7000 km 


a = 7600 km 


a = 7000 km 


Sfi (EIGEN-GRACE02S-GGM02S) 


14% 


20% 


4% 


8% 


11% 


17% 


5fi (EIGEN-GRACE02S-GGM03S) 


32% 


45% 


10% 


18% 


26% 


40% 


5fi (EIGEN-GRACE02S-ITG-Grace02s) 


33% 


43% 


9% 


17% 


26% 


38% 


5fi (EIGEN-GRACE02S-ITG-Grace03s) 


31% 


44% 


9% 


17% 


26% 


40% 


Sn (GGM02S-GGM03S) 


22% 


31% 


7% 


13% 


18% 


27% 


5fi (GGM02S-ITG-Grace02s) 


22% 


32% 


7% 


13% 


17% 


26% 


5fi (GGM02S-ITG-Grace03s) 


21% 


32% 


7% 


13% 


17% 


27% 


SfJi (GGM03S-ITG-Grace02s) 


5% 


7% 


2% 


3% 


5% 


7% 


Sfi (GGM03S-ITG-Grace03s) 


6% 


8% 


2% 


3% 


5% 


7% 


5n (ITG-Grace03s-ITG-Grace02s) 


7% 


10% 


3% 


5% 


7% 


10% 



Table 3: Systematic percent error Sfi in the measurement of the Lense-Thirring effect with the combination of 
eq. ([2]) according to eq. (|8|) and the difference AC^o among the best estimates for the even zonal coefficients for 
several Earth's gravity solutions up to degree t = 20 for different values of the inclination i of LARES close to 
polar geometry and for a = 7600 km, a = 7000 km (e = 0.001). See also Figure [3] and Figured! 





i = 85 deg 


i = 8 


9 deg 


i = 91 deg 


a = 7600 km 


a = 7000 km 


a = 7600 km 


a = 7000 km 


a = 7600 km 


a = 7000 km 


Sfi (EIGEN-GRACE02S-GGM02S) 


14% 


27% 


7% 


15% 


16% 


44% 


5fi (EIGEN-GRACE02S-GGM03S) 


32% 


52% 


16% 


28% 


34% 


82% 


5fi (EIGEN-GRACE02S-ITG-Grace02s) 


30% 


50% 


15% 


28% 


32% 


77% 


5fi (EIGEN-GRACE02S-ITG-Grace03s) 


31% 


53% 


16% 


30% 


34% 


83% 


Sn (GGM02S-GGM03S) 


22% 


39% 


11% 


21% 


24% 


63% 


5fi (GGM02S-ITG-Grace02s) 


21% 


38% 


10% 


21% 


23% 


61% 


5fi (GGM02S-ITG-Grace03s) 


22% 


40% 


11% 


22% 


24% 


64% 


SfJi (GGM03S-ITG-Grace02s) 


6% 


10% 


3% 


5% 


6% 


16% 


Sfi (GGM03S-ITG-Grace03s) 


5% 


9% 


3% 


5% 


6% 


15% 


5n (ITG-Grace03s-ITG-Grace02s) 


8% 


14% 


4% 


7% 


9% 


21% 



Table 4: Systematic percent error 5/j, in the measurement of the Lense- 
Thirring effect with LAGEOS, LAGEOS II and LARES according to eq. © 
and SJ£ = AJg up to degree t = 60 for the global Earth's gravity solutions 
considered here for a = 7600, km i = 71 deg, e = 0.001. 



Models compared 


Sfi(%) 


EIGEN-GRACE02S-GGM02S 


36% 


EIGEN-GRACE02S-GGM03S 


52% 


EIGEN-GRACE02S-ITG-Grace02s 


54% 


EIGEN-GRACE02S— ITG-Grace03s 


53% 


GGM02S-GGM03S 


24% 


GGM02S-ITG-Grace02s 


28% 


GGM02S-ITG-Grace03s 


26% 


GGM03S-ITG-Grace02s 


5% 


GGM03S-ITG-Grace03s 


5% 


ITG-Grace02s-ITG-Grace03s 


4% 
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